arXiv:1506.06864vl [hep-ph] 23 Jun 2015 


ALBERTA-THY-12-15, DESY-15-100, SFB/CPP-14-127, TTP15-021, TUM-HEP-1001/15 


Next-to-next-to-next-to-leading order QCD prediction for the top anti-top S-wave 
pair production cross section near threshold in e+e“ annihilation 

Martin Beneke/ Yuichiro Kiyo,^ Peter Marquard,^ Alexander Penin,^’^ Jan Piclum,® and Matthias Steinhauser® 

^ Physik Department T31, James-Franck-Strafie 1, 

Technische Universitdt Miinchen, D-8574-8 Garching, Germany 
‘^Department of Physics, Juntendo University, Inzai, Ghiba 270-1695, Japan 
^Deutsches Elektronen Synchrotron DESY, Platanenallee 6, D-15738 Zeuthen, Germany 
'^Department of Physics, University of Alberta, Edmonton AB TOG 2J1, Canada 
^ Institut fur Theoretische Teilchenphysik, Karlsruhe Institute of Teehnology (KIT), 76128 Karlsruhe, Germany 

^Albert Einstein Center for Fundamental Physics, 

Institute for Theoretical Physics, Sidlerstrasse 5, C 11-3012 Bern, Switzerland 

(Dated: June 17, 2015) 

We present the third-order QCD prediction for the production of top-anti-top quark pairs in electron- 
positron collisions close to the threshold in the dominant S-wave state. We observe a significant 
reduction of the theoretical uncertainty and discuss the sensitivity to the top quark mass and width. 
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Among the main motivations for building a future high- 
energy electron-positron collider in steps of increasing 
center-of-mass energy are precise measurements at -y/s Ri 
345 GeV close to the production threshold of top anti¬ 
top quark pairs. The peculiar behaviour of the cross sec¬ 
tion allows for the precise determination of a number of 
Standard Model parameters, most prominently the top 
quark mass. To date the most precise measurement of 
nit = 173.34±0.27(stat)±0.71(syst) GeV comes from the 
hadron colliders Fermilab Tevatron and CERN LHG [I| 
and is based on the reconstruction of the top and anti¬ 
top quarks through their decay products. This approach 
and the value quoted above are plagued by unknown rela¬ 
tions of the extracted mass value m* to top quark masses 
in the pole or MS renormalization scheme, which may 
well exceed 1 GeV. At hadron colliders there are also 
methods to determine directly a well-defined top quark 
mass, such as the extraction of mt from top quark cross 
section measurements. However, the final precision is of 
the order of a few GeV and thus significantly worse. At 
an electron-positron collider, on the other hand, scans 
of the top anti-top pair production threshold can lead 
to very precise measurements of well-defined mass val¬ 
ues with a statistical accuracy of only 20-30 MeV ii- 
Besides the top quark mass also its decay width and the 
strong coupling constant can be extracted with an accu¬ 
racy of 21 MeV 0 and 0.0009 0, respectively. A recent 
study has shown that for a Higgs boson with a mass of 
about 125 GeV the top quark Yukawa coupling can be 
obtained with a statistical uncertainty of only 4.2% Q. 
These numbers pose several challenges to theory. 

A crucial input to reach the aimed precision is a pre¬ 
cise calculation of the top anti-top pair production cross 
section in the threshold region. While the fundamen¬ 
tal theory of quantum chromodyanmics (QCD) is well- 
established, performing calculations of quantum correc¬ 
tions to the very high accuracies demanded here is very 
difficult indeed. The problem is further complicated by 
the fact that in the threshold region the colour Coulomb 


potential oc ag/r, where as denotes the strong coupling, 
can no longer be treated as a perturbation even though 
CTs “C 1- Standard perturbation theory in as breaks down 
and resummation is required. 

The relevant techniques have been developed in the 
1990s in the framework of effective field theory (EFT), 
which accounts for the different dynamical scales in the 
problem. For a heavy quark anti-quark system at thresh¬ 
old there are three relevant scales, the hard scale m, the 
potential and soft scale mv, and the ultrasoft scale mv'^, 
where m denotes the mass of the quark and v its velocity. 
Since r; <C 1 there is indeed a strong hierarchy between 
these scales and thus it is possible to construct a tower of 
effective theories taking QCD as starting point. In a first 
step one arrives at non-relativistic QCD (NRQCD) d, Q 
by integrating out the hard modes. Afterwards all other 
modes are integrated out except the ones which are 
needed to describe a physical non-relativistic quark anti¬ 
quark system—potential modes for the quarks with en¬ 
ergy 0{mv^) and three-momentum 0{mv), and ultrasoft 
gluons with four-momentum 0{mv‘^). The corresponding 
theory is potential NRQCD (PNRQCD) d]. In practice, 
the separation of the different modes and scales is done 
with the threshold expansion of Feynman diagrams [ 3 ]. 
The concepts and tools required to perform the computa¬ 
tion of the cross section near threshold with the accuracy 
reported in this Letter are summarized in Q. 

Within the EFT the normalized total cross section can 
be written in the form (see, e.g., @) 
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where Cy and dy denote NRQCD matching coefficients, 
E = y/s — 2mt, G{E) is the non-relativistic two-point 
Green function, and (Tq = 47ra^/(3s) is the cross section 
for the production of a fJ-~ pair in the limit of large 
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center-of-mass energy ^/s. The ellipsis refers to terms 
which are beyond the next-to-next-to-next-to-leading or¬ 
der (N^LO) in the expansion in as and v. 

The next-to-next-to-leading order (NNLO]_QCD cor¬ 
rections were computed in the late 1990s the 

results of several groups were sumirrarized and compared 
in [^. The first NNLO calculations were expressed in 
terms of the top quark pole mass and found large correc¬ 
tions to the location of the cross section peak near thresh¬ 
old casting doubt on the possibility to perform a very 
accurate mass measurement. It was pointed out in [l6| 
that these corrections are an artifact of the renormal¬ 
ization convention, which could be avoided by choosing a 
scheme that is less sensitive to uncalculable long-distance 
effects of QCD. Subsequent calculations of the top anti¬ 
top threshold all employ mass renormalization conven¬ 
tions different from the on-shell scheme. Here we use the 
potential-subtracted (PS) mass [l3]. However, even with 
this improvement, large corrections to the height of the 
cross sectioir peak are observed at NNLO, which moti¬ 
vates the N^LO QCD calculation, the result of which is 
reported in this Letter. 

There has been quite some effort to resum logarith¬ 
mic terms in the velocity of the produced top quarks 
and obtain so-called next-to-next-to-leading logarithmic 
(NNLL) approximations [H, [ip. The most complete 
analysis has been performed in [20|, where new ultrasoft 
terms have been included. The NNLL approximation 
already contains the Inu enhanced terms of the N^LO 
correction, but not the “constant” terms. Partial results 
have shown that these constant terms are as large as the 
logarithmic terms, at least in individual pieces of the cal¬ 
culation [2ll - l^ . 

Since more than ten years several groups have com¬ 
puted building blocks of the N^LO correction, which 
can be subdivided into matching coefficients of NRQCD 
and PNRQCD and higher order corrections to G{E). 
In Q a detailed discussion of the individual contribu¬ 
tions can be found. Among them are the three-loop 
corrections to the static potential (23 - 1^ . certain 0{e) 
terms of hi^er-order potentials in dimensional regular¬ 
ization JS [2^ [2^, and three-loop fermionic corrections 
to Cy 1291 l30|. Furthermore, ultrasoft corrections to 
G{E) 1^ and all Coulombic contributions up to the third 
order |^. Recently, the last two building blocks have 
been computed, which allows us to put together for the 
first time a prediction for the cross section near threshold 
with N^LO accuracy. First, the purely gluonic three-loop 
correction to the QCD-to-NRQCD vector current match¬ 
ing coefficient Cy appearing in 

+ • ■ •, ( 2 ) 

QlTlt 

has been computed [l^, completing the matching calcu¬ 
lations up to the so-called sin glet diagrams. (In other 
contexts and at second order these singlet contri¬ 
butions have been shown to be small.) Second, the 
calculation of the PNRQCD two-point function G{E) 


of currents been completed to third order 

in PNRQCD perturbation theory with the computation 
of the single and double non-Coulomb potential inser¬ 
tions [ 2 ^. We note that PNRQCD includes the leading- 
order Coulomb potential in the unperturbed Lagrangian. 
Perturbation theory around this interacting Lagrangian 
amounts to the resummation of the standard loop expan¬ 
sion in as- The various contributions have been encoded 
in a Mathematica program. 

There is a number of non-QCD effects, which are ex¬ 
pected to be relevant to a realistic cross section predic¬ 
tion. a) The existence of a Higgs boson affects the pro¬ 
duction process and generates an additional short-range 
potential, b) The same holds for electroweak and electro¬ 
magnetic corrections, amongst which is the electrostatic 
potential between the charged top quarks, c) The rapid 
electroweak decay t —>■ W^b of the top quark is respon¬ 
sible for the fact that the distinct toponium bound-state 
poles are smeared to a broad resonance near threshold 
(see Fig.[T]below). For the same reason, it will be manda¬ 
tory to consider the process e+e“ —>■ W^W~bb includ¬ 
ing the non-resonant contributions. Non-resonant effects 
have been computed to NLO and partially to NNLO ac¬ 
curacy [ssl - l^ . and can be naturally separated in the 
framework of unstable particle EFT [13, d) Elec¬ 
tromagnetic initial-state radiation generates large loga¬ 
rithms ln(4m(/mg), where me is the electron mass, which 
must be summed, e) Eurthermore, there are collider- 
specific effects such as beamstrahlung and the luminosity 
spectrum Q. Neither of these effects is discussed further 
in the following. The purpose of the present Letter is to 
demonstrate that the third-order QCD calculation brings 
the QCD contributions under control, which is a prereq¬ 
uisite for all further studies, including the non-QCD ef¬ 
fects above. We refer to [1^ for a study of the Higgs con¬ 
tribution and the prospects for determining the Yukawa 
coupling, and the impact of the NLO non-resonant cor¬ 
rection on the top mass determination. 

We now turn to the discussion of the N^LO QCD re¬ 
sult. In order to better compare the successive orders we 
do not include here the known contribution [^ from the 
axial-vector Z-boson coupling in the production process 
e+e” —>■ Z* —?► tt, which starts at NNLO and is below 
1%. Our result therefore refers to the S-wave production 
cross section. Unless stated otherwise we use the follow¬ 
ing input values for the top quark mass in the PS scheme, 
the top quark width and the strong coupling constant: 

= 20GeV) = I7I.5 GeV, F; = 1.33 GeV, 
a^(Mz) = 0.1185 ±0.0006. (3) 

Furthermore, the renormalization scale is varied between 
50 and 350 GeV. Note that an unstable behaviour of 
the perturbative series is expected for scales below fi Ri 
40 GeV m, [sij. The scale related to the separa¬ 
tion of resonant and non-resonant contributions to the 
e+e“ —>■ W^W~bb process is fixed to ^y, = 350 GeV. 
The dependence on this choice is numerically negligible. 

The main result of this Letter is shown in Fig.[l] where 
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FIG. 1. Scale dependence of the cross section near thresh¬ 
old. The NLO, NNLO and N®LO result is shown in blue, red 
and black, respectively. The renormalization scale is varied 
between 50 and 350 GeV. 


the total cross section is shown as a function of the center- 
of-mass energy ^/s. The previous NLO and NNLO pre¬ 
dictions are also shown for comparison to the new N^LO 
result (black, solid). The bands are obtained by varia¬ 
tion of the renormalization scale in the specified range. 
After the inclusion of the third-order corrections one ob¬ 
serves a dramatic stabilization of the perturbative predic¬ 
tion, in particular in and below the peak region. In fact, 
the N^LO curve is entirely contained within the NNLO 
one. This is different above the peak position where a 
clear negative correction is observed when going from 
NNLO to N^LO. For example, 3 GeV above the peak 
this amounts to —8%. This arises from the large negative 
three-loop correction to the matching coefficient Cy . 

The theoretical precision of the third-order QCD result 


as measured by the residual scale dependence is high¬ 
lighted in Fig. [5J which shows normalized to a ref¬ 
erence prediction defined at /r = 80 GeV. The width of 
the shaded band corresponds to an uncertainty of about 
±3% with some dependence on the center-of-mass energy 
i/s. The figure also shows the sensitivity to the top-quark 
width. The two solid lines refer to the cross section with 
Ft changed by ±100 MeV to 1.43 and 1.23 GeV, respec¬ 
tively, computed with /i = 80 GeV and normalized to 
the reference prediction. Decreasing the width implies 
a sharper peak, i.e. an enhancement in the peak region, 
and a suppression towards the non-resonant region below 
the peak. A few GeV above the peak the cross section 
is largely insensitive to the width. Increasing the width 
leads to the opposite effects. This pattern is clearly seen 
in Fig. [2J which also demonstrates that a ±100 MeV de¬ 
viation from the width predicted in the Standard Model 
leads to a cross section change near and below the peak 
far larger than the uncertainty from scale variation. 

We now turn to the question to what accuracy the 
top quark mass can be determined. Even if we focus 
only on the theoretical accuracy, a rigorous analysis re¬ 
quires accounting for the specifics of the energy points 
of the threshold scan and the correlations. However, a 
good indication is already provided by looking at the po¬ 
sition and height of the resonance peak. Fig.[3]shows this 
information at LO, NLO, NNLO and N^LO, where the 
outer error bar reflects the uncertainty due to the renor¬ 
malization scale and as variation, added in quadrature, 
and the inner error bar only takes the Og uncertainty 
into account. The central point refers to the value at 
the reference scale fi = 80 GeV. There is a relatively big 
jump from LO to NLO of about 310 MeV, approximately 
150 MeV from NLO to NNLO, which reduces to only 
64 MeV from NNLO to N^LO. Furthermore, the NNLO 
and N^LO uncertainty bars show a significant overlap. 
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FIG. 2. Scale dependence (hatched area) of the N^LO cross 
section relative to the reference prediction. Overlaid are pre¬ 
dictions for two different values of Ft, again normalized to the 
reference prediction. See text for details. 


FIG. 3. Position and height of the cross section peak at LO, 
NLO, NNLO and N^LO. The unbounded range of the LO 
error bars to the right and up are due to the fact that the 
peak disappears for large values of the renormalization scale. 


















4 



(GeV) 

FIG. 4. Relative cross section variation when the top quark 
mass is changed by ±50 and ±100 MeV, superimposed on the 
scale dependence of the N^LO cross section prediction. See 
text for details. 


Taking into account only the uncertainty from scale vari¬ 
ation the uncertainty of the peak position amounts to 
±60 MeV at N^LO with a factor two improvement rel¬ 
ative to NNLO. The improvement is even larger for the 
peak height, which is relevant to the top quark width and 
Yukawa coupling determination as discussed above and 
in . Note that these conclusions refer to the top quark 
PS mass (and not the pole mass), and correspondingly to 
the MS mass, which can be related to the PS mass with 
an accuracy of about 20 MeV (4ll |. 

We display the sensitivity of the cross section to the 
top quark mass in Fig. SI which is the same as Fig. SJ 
except that the curves superimposed to the relative scale 
variation of the N^LO cross section now refer to shifts 
of the top PS mass by ±50 and ±100 MeV. The shape 
of these curves is easily understood, as a change of the 
mass by 5mt mainly shifts the peak of the cross section 
and is largely equivalent to a shift of ^ by 26mt- The 
figure demonstrates that the largest sensitivity to the 
mass occurs around 1.5 GeV below the peak. A varia¬ 
tion of the mass by ±50 MeV changes the cross section 


at -v/i = 342.5 GeV by ±6%, compared to a scale uncer¬ 
tainty of ±3.8%. Accounting for the characteristic shape 
of variation, we may conclude that the theoretical un¬ 
certainty on the top quark mass should be well below 
±50 MeV. 

Comparing the N^LO result to the NNLO results with 
Inu resummation BUI, we notice that in both ap¬ 
proaches the correction to the peak height relative to 
NNLO is negative (with the default scale choice) [l^ . 
while the theoretical uncertainty of the cross section nor¬ 
malization is reduced to about ±3% compared to the 
±5% quoted in [^ . 

In summary, we presented the third-order QCD calcu¬ 
lation of the top quark production cross section in e+e" 
annihilation in the threshold region. Despite the extra 
complication of non-relativistic resumination to all or¬ 
ders in QCD perturbation theory, this is one of only a 
few collider processes now known at N^LO in QCD. The 
third-order calculation leads to a large reduction of the 
QCD theoretical uncertainty to about ±3%, and thereby 
solves a long-standing issue regarding the reliability of 
the QCD prediction for this process. Top quark mass de¬ 
terminations with theoretical errors below 50 MeV now 
appear feasible. Further improvement by going to the 
next order in QCD is currently unrealistic. However, it 
would be desirable to include the information about log¬ 
arithmic effects in v beyond N^LO already contained in 
the NNLL computations [illig. More importantly, with 
QCD effects under control as demonstrated here, further 
studies of non-QCD effects are now well motivated—and 
required. 

ACKNOWLEDGMENTS 

This work has been supported by DFG Sonderforschungs- 
bereich Transregio 9 “Coinputergestiitzte Theoretische 
Teilchenphysik”, and the DFG Gottfried Wilhelm Leib¬ 
niz programme. The work of Y.K. was supported in part 
by Grant-in-Aid for Scientific Research No. 26400255 
from MEXT, Japan. P.M. was supported in part by the 
European Commission through contract PITN-GA-2012- 
316704 (HIGGSTOOLS). The work of A.P. is partially 
supported by NSERC. 


[1] ATLAS and CDF and CMS and DO Collaborations, 
arXiv:1403.4427 [hep-ex]. 

[2] K. Seidel, F. Simon, M. Tesar and S. Poss, Eur. Phys. J. 
C 73, 2530 (2013), arXiv:1303.3758 [hep-ex]. 

[3] T. Horiguchi, A. Ishikawa, T. Suehara, K. Fujii, Y. Sum- 
ino, Y. Kiyo and H. Yamamoto, arXiv:1310.0563 [hep- 
ex]. 

[4] B. A. Thacker and G. P. Lepage, Phys. Rev. D 43 (1991) 
196. 

[5] G. P. Lepage, L. Magnea, C. Nakhleh, U. Magnea 
and K. Hornbostel, Phys. Rev. D 46 (1992) 4052 


hep-lat/9205007 . 

[6] A. Pineda and J. Soto, Nucl. Phys. Proc. Suppl. 64, 428 
(1998) hep-ph/9707481 . 

[7] M. Beneke and V. A. Smirnov, Nucl. Phys. B 522, 321 
(1998) hep-ph/9711391 . 

[8] M. Beneke, Y. Kiyo and K. Schuller, arXiv: 1312.4791 
[hep-ph]. 

[9] A. H. Hoang and T. Teubner, Phys. Rev. D 58, 114023 
(1998) hep-ph/9801397 . 

[10] K. Melnikov and A. Yelkhovsky, Nucl. Phys. B 528, 59 
(1998) hep-ph/9802379 . 








5 


[11] A. A. Penin and A. A. Pivovarov, Nucl. Phys. B 550, 
375 (1999) hep-ph/9810496 . 

[12] M. Beneke, A. Signer and V. A. Smirnov, Phys. Lett. B 
454, 137 (1999) hep-ph/9903260 . 

[13] A. H. Hoang and T. Teubner, Phys. Rev. D 60, 114027 
(1999) hep-ph/9904468 . 

[14] T. Nagano, A. Ota and Y. Sumino, Phys. Rev. D 60, 
114014 (1999) hep-ph/9903498 . 

[15] A. H. Hoang, M. Beneke, K. Melnikov, T. Nagano, 
A. Ota, A. A. Penin, A. A. Pivovarov and A. Signer et 
ai, Eur. Phys. J. direct C 2, 1 (2000) hep-ph/0001286 . 

[16] M. Beneke, in; Proceedings of the 33rd Rencontres de 
Moriond: Electroweak Interactions and Unified Theories, 
Les Arcs, France, 14-21 Mar 1998, hep-ph/9806429 

[17] M. Beneke, Phys. Lett. B 434, 115 (1998) 

hep-ph/9804241 . 

[18] A. H. Hoang, A. V. Manohar, 1. W. Stewart and T. Teub¬ 
ner, Phys. Rev. D 65, 014014 (2002) hep-ph/0107144 . 

]19] A. Pineda and A. Signer, Nucl. Phys. B 762, 67 (2007) 
hep-ph/0607239 . 

[20] A. H. Hoang and M. Stahlhofen, JHEP 1405 (2014) 121, 
arXiv:1309.6323 [hep-ph]. 

[21] M. Beneke and Y. Kiyo, Phys. Lett. B 668 (2008) 143, 
arXiv:0804.4004 [hep-ph]. 

[22] P. Marquard, J. H. Piclum, D. Seidel and M. Steinhauser, 
Phys. Rev. D 89 (2014) 3, 034027, arXiv:1401.3004 [hep- 
ph]. 

[23] M. Beneke, Y. Kiyo and K. Schuller, “Third-order cor¬ 
rection to top-quark pair production near threshold H. 
Potential contributions”, to be published. 

[24] A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. 
Lett. B 668, 293 (2008), arXiv:0809.1927 [hep-ph]. 

[25] A. V. Smirnov, V. A. Smirnov and M. Steinhauser, Phys. 
Rev. Lett. 104, 112002 (2010), arXiv:0911.4742 [hep-ph]. 

[26] C. Anzai, Y. Kiyo and Y. Sumino, Phys. Rev. Lett. 104, 
112003 (2010), arXiv:0911.4335 [hep-ph]. 


[27] S. Wiister, Heavy quark potential at order a^fm^ [in 
German], Diploma Thesis RWTH Aachen University 
(2003). 

[28] M. Beneke, Y. Kiyo, P. Marquard, A. Penin, J. Pi¬ 
clum, D. Seidel and M. Steinhauser, Phys. Rev. Lett. 
112 (2014) 15, 151801, arXiv:1401.3005 [hep-ph]. 

[29] P. Marquard, J. H. Piclum, D. Seidel and M. Steinhauser, 
Phys. Lett. B 678, 269 (2009), arXiv:0904.0920 [hep-ph]. 

[30] P. Marquard, J. H. Piclum, D. Seidel and M. Steinhauser, 
Nucl. Phys. B 758, 144 (2006) hep-ph/0607168 . 

[31] M. Beneke, Y. Kiyo and K. Schuller, Nucl. Phys. B 714, 
67 (2005) hep-ph/0501289 . 

[32] B. A. Kniehl, A. Onishchenko, J. H. Piclum and M. Stein¬ 
hauser, Phys. Lett. B 638, 209 (2006) hep-ph/0604072 . 

[33] M. Beneke, B. Jantzen and P. Ruiz-Femenia, Nucl. Phys. 
B 840, 186 (2010), arXiv:1004.2188 [hep-ph]. 

[34] A. A. Penin and J. H. Piclum, JHEP 1201 (2012) 034, 
arXiv:1110.1970 [hep-ph]. 

[35] B. Jantzen and P. Ruiz-Femenia, Phys. Rev. D 88, 
054011 (2013), arXiv: 1307.4337 [hep-ph]. 

[36] P. Ruiz-Femenia, Phys. Rev. D 89 (2014) 9, 097501, 
arXiv:1402.1123 [hep-ph]. 

[37] M. Beneke, A. P. Chapovsky, A. Signer and 

G. Zanderighi, Phys. Rev. Lett. 93 (2004) 011602 
hep-ph/0312331 . 

[38] M. Beneke, A. P. Ghapovsky, A. Signer and G. Zan¬ 
derighi, Nucl. Phys. B 686 (2004) 205 hep-ph/0401002 . 

[39] M. Beneke, A. Maier, J. Piclum and T. Rauh, to be pub¬ 
lished. 

[40] M. Beneke, J. Piclum and T. Rauh, Nucl. Phys. B 880 
(2014) 414, arXiv:1312.4792 [hep-ph]. 

[41] P. Marquard, A. V. Smirnov, V. A. Smirnov and 
M. Steinhauser, Phys. Rev. Lett. 114 (2015) 14, 142002, 
arXiv: 1502.01030 [hep-ph]. 






